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The two-dimensional freezing transition is very different from its three-dimensional counterpart. 
In contrast, the glass transition is usually assumed to have similar characteristics in two and three 
dimensions. Using computer simulations we show that glassy dynamics in supercooled two- and 
three-dimensional fluids are fundamentally different. Specifically, transient localization of particles 
upon approaching the glass transition is absent in two dimensions, whereas it is very pronounced 
in three dimensions. Moreover, the temperature dependence of the relaxation time of orientational 
correlations is decoupled from that of the translational relaxation time in two dimensions but not in 
three dimensions. Lastly, the relationships between the characteristic size of dynamically heteroge¬ 
neous regions and the relaxation time are very different in two and three dimensions. These results 
strongly suggest that the glass transition in two dimensions is different than in three dimensions. 


INTRODUCTION 

In two-dimensional (2D) solids, thermal fluctuations 
destroy crystalline order, displacement correlations in¬ 
crease logarithmically, and density correlations decay ac¬ 
cording to power laws mm- However, there can be long- 
range bond-orientational order in 2D [2]. The transition 
from the 2D fluid phase to the solid phase can occur in 
two steps with an intermediate phase characterized by 
an exponential decay of the density correlations and a 
power-law decay of the bond-orientational correlations 
mm- In contrast, in three-dimensional (3D) solids fluc¬ 
tuations do not destroy crystalline order [3], and long- 
range translational and rotational order emerge together 
at the freezing transition. 

Despite these differences between two- and three- 
dimensional ordered solids, the formation of an amor¬ 
phous solid upon supercooling a fluid, i.e. the glass tran¬ 
sition, is generally assumed to have similar characteris¬ 
tics in 2D and 3D . This assumption is reflected in the 
trivial dimensional dependence of most glass transition 
theories [B]. 

We show that structural relaxation of supercooled flu¬ 
ids in two dimensions is different than in three dimen¬ 
sions. While we find the transient localization often as¬ 
sociated with glassy dynamics in three dimensions, we do 
not find any transient localization in two dimensions if 
we simulated systems large enough to remove any finite 
size effects. Furthermore, the temperature dependence 
of the bond-orientational correlation time is decoupled 
from that of the translational relaxation time time in two 
dimensions, but these relaxation times have very simi¬ 
lar temperature dependence in three dimensions. Along 
with these differences in structural relaxation, we also 
find that the characteristic size of regions of correlated 
mobility, dynamic heterogeneities, increases faster with 
the structural relaxation time in two dimensions than in 
three dimensions, and these regions are more ramified in 
two dimensions than three dimensions. Lastly, we show 
that the structural relaxation and heterogeneous dynam¬ 


ics depends on the underlying dynamics in two dimen¬ 
sions. 


RESULTS 

Structural relaxation. To demonstrate the differ¬ 
ences between glassy dynamics in two and three dimen¬ 
sions, we focus on two closely related glass-forming fluids: 
the 3D 80:20 binary Lennard-Jones system introduced by 
Kob and Andersen [7], and its 2D variant that has the 
same interaction potentials but a 65:35 composition to 
avoid crystallization [8]. To simulate the relaxation in 
these systems we used the standard Newtonian dynam¬ 
ics [5]. We also simulated the 2D system using Brownian 
dynamics and we comment on the differences between 
these two dynamics. See Methods for the simulation de¬ 
tails and the reduced units used to present the results. 
We examined three other 2D glass formers and one ad¬ 
ditional 3D glass former (see Methods for details of the 
systems) and present results, which are qualitatively the 
same as the results for the KA system, in the supplemen¬ 
tal material. 

In 3D, the dominant feature in the dynamics of deeply 
supercooled glass-forming fluids is transient localization 
of individual particles [B] , which is illustrated in the inset 
to Fig. The transient localization results in character¬ 
istic plateaus of the self-intermediate scattering function, 
Fsik-,t) = AT-i (k=7.2 in 3D and 

6.28 in 2D), shown in Fig.f^a), and the mean-square dis¬ 
placement, (Sr'^it)) = N-^{J2n[TCn(t) - r„(0)]^), shown 
in Fig. The plateaus extend to longer and longer 
times upon approaching the glass transition. Similar 
plateaus are observed in the collective scattering func¬ 
tion, F{k-,t) = (^X^n m '"^hich de¬ 

scribes relaxation of the density field (not shown), and in 
the correlation function quantifying bond-orientational 
correlations in 3D, CQ{t), shown in Fig. B: (see Methods 
for the definition of Cq). Qualitatively similar slowing 
down of the translational and bond-orientational relax- 
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FIG. 1: Structural relaxation in two and three dimensions, a, The self-intermediate scattering function Fs{k-,t) for 
the 3D glass-former for T = 1.0, 0.8, 0.7, 0.6, 0.55, 0.5, and 0.45 listed from left to right. The mode-coupling temperature 
~ 0.435. The inset shows a trajectory plot of one small particle at T = 0.45 where the color of the trajectory plot changes 
when the particle moves more than one large particle diameter a (length of the black line) over a time of O.Itq. b. The mean 
square displacement for the 3D system showing the same temperatures as in (a), c, The bond-orientational correlation 

function Gqil) for the 3D glass-former showing the same temperatures as in (a), d. The self-intermediate scattering function 
Fs{k\t) for the 2D glass-former for T = 1.0, 0.8, 0.7, 0.6, 0.5, and 0.45 listed from left to right. The inset shows a trajectory 
plot of a small particle, and no sudden jumps are observed, e. The mean square displacement ((5r^(t)) for the 2D system 
showing the same temperatures as in (d). f. The bond-orientational correlation function for the 2D glass-former showing 

the same temperatures as in (d). 



FIG. 2: Bond angle and translational relaxation times, a, The ratio of the relaxation time for the bond-orientational 
correlation functions re and the self-intermediate scattering function Ta for the 2D system (circles) and the 3D system (squares). 
The ratio tbItc, for the 3D system is approximately constant and equal to 0.1-0.2 over the entire range of temperatures. The 
green circles are for Newtonian dynamics (ND) and the blue circles are results for Brownian dynamics (BD). b. The self¬ 
intermediate scattering function Fs{k;t) (red solid lines) and the bond-angle time correlation function Gqil) (blue lines) 
rescaled by for the 3D system, c. The self-intermediate scattering function Fs{k\t) (green solid lines) and the bond-angle 
time correlation function G\sj{t) (violet dashed lines) rescaled by Tc for the 2D system. 


ation in 3D glass forming fluids is analogous to the simul¬ 
taneous appearance of translational and rotational long 
range order in 3D crystalline solids. 

The transient localization observed in 3D glassy dy¬ 
namics is absent in 2D, as showed in the inset to Fig.[^. 


Correspondingly, there is no intermediate time plateau in 
the self-intermediate scattering function in 2D, Fig. [^. 
The final decay of F's(fc; t), which in 3D is well described 
by a stretched exponential, is replaced by a very slow 
decay in 2D. The intermediate time plateau in the mean- 
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square displacement observed in 3D is replaced by an 
extended sub-diffusive regime in 2D, Fig. [^. However, 
an intermediate time plateau is observed in the correla¬ 
tion function quantifying bond-orientational correlations 
in 2D, Fig. (see Methods for the definition of 

C^). Shown in Supplemental Figure la-c are Fs{k;t) 
and Cqi{t) for three additional 2D glass formers and 
they behave similarly. Qualitatively different behavior 
of the translational and bond-orientational relaxation in 
2D glass forming fluids is analogous to the absence of the 
translational and the presence of the bond-orientational 
long-range order in 2D solids [T]. 

To quantify decoupling between translational and 
bond-orientational relaxation we compare the tempera¬ 
ture dependence of the relaxation times characterizing 
Fs{k;t) and C(Q_^)(t), where Q and refer to 3D and 
2D correlation functions. We define the translational re¬ 
laxation time Ta through the relation Fs{k;Ta) = e~^ 
and the bond-orientational relaxation time Tg through 
C{Q^^){Tg) = e~^. At the highest temperatures the ra¬ 
tio TgjTa is less than one for both the 3D and the 2D 
glass-former, Fig. However, this ratio stays approxi¬ 
mately constant with decreasing temperature for the 3D 
glass-former, but grows monotonically for the 2D glass- 
former. In Supplemental Figure Id we show this ratio for 
the other 2D glass formers and show that the decoupling 
is a general feature of 2D glassy dynamics. In addition, in 
Figs. [^,c we show that the final translational and orien¬ 
tational relaxation satisfies the time-temperature super¬ 
position in 3D but not in 2D, and we show corresponding 
figures for an additional glass former in 2D and 3D in 
Supplemental Figure 2. Fig. clearly demonstrates the 
decoupling of the temperature dependence of the trans¬ 
lational and bond-orientational relaxation times in 2D. 

Dynamic heterogeneities. The non-exponential de¬ 
cay of Fs{k;t) is frequently attributed to the emergence 
of domains, referred to as dynamic heterogeneities, in 
which the relaxation is spatially correlated and signifi¬ 
cantly different (faster or slower) than the average relax¬ 
ation. While we find non-exponential decay in Fs(k] t) for 
3D and 2D glass-formers, the nature of the decay is very 
different and this difference is mirrored by differences in 
the heterogeneous dynamics. 

Shown in Figs. |^-d are displacement maps showing 
the center of a four million particle simulation in 2D at 
T = 0.45. The maps are created by coloring the parti¬ 
cles, whose position is shown on t = 0, according to the 
magnitude of their displacements |r„(f) — r„(0)| at a time 
t. The red particles have moved a distance equal to or 
greater than the diameter of a larger particle. There are 
large domains of particles that have moved less than a 
particle diameter even at t = 10 000. 

Considering the large dynamically heterogeneous re¬ 
gions in Figs, l^-d, it is unsurprising that we also find 
large finite size effects. Shown in Fig. is Fs{k;t) cal¬ 
culated for different size systems at the same temper- 
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FIG. 3: Dynamic heterogeneities, a-d. Displacement 
maps of the center of a system of 4 million particles that 
show the position of the particles at t = 0 colored by the 
displacement of the particle at a later time t for T = 0.45. 
There are approximately 250 000 particles in each map. The 
scale bar in d corresponds to 100 larger particle diameters, e, 
The self-intermediate scattering function Fs{k-,t) calculated 
for systems of 100 to 4 million particles. There are clear fi¬ 
nite size effects for less than one million particles. 


ature as shown in Figs. -d. A plateau reminiscent 
of the plateau in 3D systems is present for the smaller 
systems but gradually disappears with increasing system 
size. Similar finite size effects are also evident in the mean 
square displacement. Supplemental Figure 3a, and the 
inherent structure dynamics. Supplemental Figure 3b. 

To quantify dynamic heterogeneity shown in 
Fig. ^ -d we use a four-point structure factor 
S 4 {q;t) [T5] constructed from overlap functions 

Wn{a;t) = 0(a - |r„(t) - r„(0)|), where ©(•) is 
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FIG. 4: Properties of dynamic heterogeneities, a, The dynamic correlation length ^ 4 ,{Ta) versus the translational relax¬ 
ation time Ta- The circles are the results for the 2D system where the underlying dynamics is Newtonian (green) and Brownian 
(blue). The lines are fits to a power law f 4 (ra) = arf. The red squares are results for the 3D system where the underlying 
dynamics is Newtonian. The red line is a fit to ^4,{Tai) = 6[ln(ra/ro)]^^®. b, The dynamic susceptibility versus the dy¬ 

namic correlation length ^ 4 (ra) for the 2D system where the underlying dynamics are Newtonian (green circles) and Brownian 
(blue circles). The red squares are results for the 3D system. The lines are fits to the power laws ~ ^ 4 {Ta)^ in 3D and 

X 4 {Ta) ~ ^ 4 ,{TaY'^ in 2D. c, The relaxation time Ta and the dynamic correlation length ^ 4 (Ta) where the underlying dynamics 
are Newtonian (green) and Brownian (blue) for the 2D system, and Newtonian for the 3D system (red). 


Heaviside’s step function. The parameter a is chosen 
such that {J2n'Wn{a',Ta)) « Fs{k;Ta), which 

results in a = 0.25 in 3D and a = 0.22 in 2D. To 
characterize the slow domains we calculate S' 4 (g;t) = 
{J2n Wnia; t)wm{a; (note 

that w„(a;t) restricts the sums over the particles that 
moved less than a over a time t). The characteristic 
size of dynamically heterogeneous regions is quantified 
through the dynamic correlation length ^ 4 ( 1 ), which 
is determined from fitting S 4 {q;t) for small q to the 
Ornstein-Zernicke form X 4 (t )/{1 + [( 7 ^ 4 (t)]^}. Here X 4 (t) 
is the dynamic susceptibility, which characterizes the 
overall strength of the dynamic heterogeneity. 

In Fig. we show the correlation between the trans¬ 
lational relaxation time, Tq,, and the dynamic correlation 
length calculated at Tq,, ^ 4 {Ta), for the 3D and 2D glass 
forming fluids. While for the 3D system we find that a 
power law is a poor description for an extended range 
of Tq, and a better description is 'C 4 ('Ti) [lii('rct/'n 3 )]^^^ 

(red line in Fig.|^), we find that a power law ^ 4 {Ta) ~ 
with P = 1.0 ± O.f describes the full range of results well 
for the 2D system. We show results for additional glass 
formers in 2D and 3D in Supplemental Figure 4. Note 
that similar power law behavior was observed in simu¬ 
lations of 2D granular fluids m- In Fig. & we show 
that the relationship between the dynamic susceptibility 
and the dynamic correlation length is fundamentally dif¬ 
ferent in 3D and 2D. For 3D systems X 4 (to,) ~ 'C 4 ('^q)^ 
at low temperatures, which implies compact dynami¬ 
cally heterogeneous regions. For 2D systems we observe 
X 4 ('ra) ^ 4 (T'ct)^’^, which suggests more ramified dynam¬ 

ically heterogeneous regions, see Fig. P. 

Dependence on the microscopic dynamics. 

Lastly, we discuss the dependence of the long-time re¬ 


laxation in 2D on the underlying microscopic dynamics 
and two important consequences. In 3D, an important 
finding is that the long-time dynamics does not depend 
on the microscopic dynamics; the same long-time dynam¬ 
ics has been observed in simulations using Newtonian [7], 
stochastic Brownian [ 12 ] and Monte Carlo [T3| dy¬ 
namics. This result can be rationalized within the mode¬ 
coupling approach mi. Surprisingly, we find that in 2D 
the long-time dynamics is quite different in the case of 
microscopic Newtonian and Brownian dynamics. The re¬ 
sults corresponding to the those shown in Fig. P -f for 
the Newtonian case are shown in Fig.j^for the Brownian 
case. Notably, the decay of Fs(fc;t) is strikingly differ¬ 
ent for the Brownian simulations than for the Newtonian 
simulations. 

Importantly, the temperature dependence of the trans¬ 
lational relaxation time is also decoupled from orienta¬ 
tional relaxation time in the case of Brownian dynamics. 
Fig. but the ratio to/to, is not as large for Brown¬ 
ian dynamics than for Newtonian dynamics. In addition, 
we find a power law relationship ^ 4 (tq,) ~ between 
the dynamic correlation length and the relaxation time 
but with P = 0.36 ± 0.05, which is a different exponent 
than obtained for Newtonian dynamics, see Fig. How¬ 
ever, we find that the relationship between the strength 
of the dynamic heterogeneity and the dynamic correla¬ 
tion length in 2D is the same for Brownian and Newto¬ 
nian dynamics, see Fig. |^. The latter two results show 
that the universality of the relationships between the re¬ 
laxation time and properties of heterogeneous dynamics 
that we found in 3D [15] is absent in 2D. Furthermore, 
a full description of heterogeneous dynamics in 2D must 
also include the influence of the microscopic dynamics, 
and descriptions solely in terms of the structure or the 
potential energy landscape are not sufficient in 2D. 
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FIG. 5: Structural relaxation in two dimensions with Brownian dynamics, a, The self-intermediate scattering function 
Fs{k-,t) for T = 1.0, 0.8, 0.7, 0.6, 0.5, 0.45 and 0.4 listed from left to right, b. The mean square displacement ((5r^(t)) showing 
the same temperatures as in (a), c. The bond angle time correlation function C'>i>(t) showing the same temperatures as in a. 
We emphasize that the 6 higher temperatures showed in this figure are the same as the temperatures showed in Figs. [i]:d)-(f). 
The hnite size effects are much more pronounced in systems evolving with Newtonian dynamics. This fact made impossible to 
simulate the 2D system at T = 0.4 with Newtonian dynamics. 


DISCUSSION 

Glassy dynamics in 2D and in 3D are profoundly dif¬ 
ferent. While we only presented detailed results for one 
glass-former, we verified that the features of the trans¬ 
lational relaxation and dynamic heterogeneity are qual¬ 
itatively the same for three additional 2D glass-formers 
(see Methods for their description and Supplemental In¬ 
formation for the results). Our results call for a re¬ 
examination of the present glass transition paradigm in 
2D. We note that there is currently no theoretical frame¬ 
work that accounts for the different dynamics observed 
in the 2D glass forming systems. However, we note that 
the dynamic picture of the Random First Order Tran¬ 
sition theory breaks down for dimensions less than two, 
and has been described as marginal for two dimensions 
[Mini- Moreover, insights gained from theoretical anal¬ 
ysis of the 2D glassy dynamics and glass transition might 
shed light onto slow dynamics and the glass transition in 
3D. It will also be interesting to investigate if the differ¬ 
ences between 2D and 3D glassy dynamics are observable 
for glass-forming fluids in confinement and at interfaces 
or surfaces, i.e. for quasi-two-dimensional systems. 

We gratefully acknowledge the support of NSF grant 
CHE 1213401. This research utilized the CSU IS- 
TeC Cray HPC System supported by NSF Grant CNS- 
0923386. 


METHODS 

Simulations 

We simulated binary mixtures of Lennard-Jones particle 
in two and three dimensions. The interaction potential is 
Vap{r) = ^tap[{(Tap/rY‘^ “ (cTc/s/r)®] where tBB = 0.5£aa, 
CAS = 1.5eAA, frss = O.SSctaa, and (jab = O.Sctaa. The re¬ 


sults are presented in reduced units where oaa = (t is the unit 
of length and caa the unit of energy. The unit of time for the 
Newtonian dynamics simulations is a’^ml^AA and the mass 
m is the same for both species. The Newtonian dynamics sim¬ 
ulations [3 were performed using LAMMPS [TS| for the 2 D 
and 3D simulations and HOOMD-blue for the 2D simulations 
m- The LAMMPS simulations were run in an NVE ensem¬ 
ble, but there is significant energy drift for the HOOMD-blue 
simulations for the lowest temperatures. Therefore, we ran 
the HOOMD-blue simulations using an NVT Nose-Hoover 
thermostat with a coupling constant r = 10. We ran at 
least one LAMMPS NVE simulation at every temperature 
to make sure that the conclusions did not depend on the 
thermostat. All the results are averages over four or more 
production runs. The equations of motion for the Brown¬ 
ian dynamics simulations [9] are r„{t) = 7 “^F„(t) -|- rjn{t), 
where 7 = 1 is the friction coefficient, F„(t) is the force 
on particle n at time t, and rjn is a random noise term. 
The random noise satisfies the fluctuation dissipation rela¬ 
tion {rin{t)rim{t')) = 2kBT^~^5{t — where 1 is the 

unit tensor. The unit of time for the Brownian dynamics 
simulation is cr^ 7 /eAA. The Brownian dynamics simulations 
were run using a modified version of LAMMPS and our in 
house developed code. 

We simulated 2D systems of 10 000 particles for T > 0.9 
and 250 000 particles for 0.5 < T < 0.8. At T — 0.45 we 
studied 4 million particles for the Newtonian dynamics simu¬ 
lations, but 250 000 particles for the Brownian dynamics simu¬ 
lations. We simulated 27 000 particles in 3D using Newtonian 
dynamics. To check that the results are independent of the 
system size for each state point for the 2D Newtonian dynamic 
simulations we ran 100 489 particle simulations and checked 
to see if the results agreed with the 250 000 particle simula¬ 
tions. At T = 0.45 they did not agree, and we increased the 
system size until we found agreement between the 4 million 
particle system and an 8 million particle system. For the 2D 
Brownian dynamics simulations we found agreement between 
10 000 particle simulations and 250 000 particle simulations 
for T > 0.45. For T = 0.4 we found that a 100 489 particle 
simulation agreed with a 250 000 particle system. 

We also examined the translational dynamics, bond- 
orientational, and dynamic heterogeneities for three addi- 
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tional systems in 2D and one additional glass forming system 
in 3D. The first system is the one studied in Ref. m and 
consists of a 32.167:67.833 binary mixture with the poten¬ 
tial Vapir) = tiaapjrY'^. The size ratios are uab = \.\uaa 
and obb = IAgaa- We simulated this system using 250 000 
particles in 2D and 100 000 particles 3D. The number den¬ 
sity p = 0.719(7^^ in 2D and p = 0.719cr^i^ in 3D. The 
second is a system introduced in Ref. [^, which consists 
of an 50:50 mixture of repulsive particles where the poten¬ 
tial Vap{r) = ^{cTapIrY'^. The size ratios are given by 
(7as = 1.2(jaa and aBB = IAoaa and the number den¬ 
sity p = N/L^ = 0.74718(t^^. We simulated 250000 par¬ 
ticles for this second additional system. The third system is 
the one introduced in Ref. [22] . which consists of an 50:50 
mixture of harmonic spheres with the interaction potential 
14/3(r) = 0.5£(1 - rjoapY for r < a^p and Vap{r) = 0 
otherwise. The size ratios are given by uab = 1-2(Jaa and 
ubb = IAoaa and p — 0.699a^\. We simulated 250 000 par¬ 
ticles for this third additional system. Some results for the 
three systems described in this paragraph are given in the 
Supplemental Material. 


Bond-Orientational Correlation Functions 

To measure bond-orientational relaxation times in 2D we 
first define where 9nm{t) is 

the angle between particle n and particle m at a time t, 
is the number of neighbors of particle n, and the sum is over 
the neighbors of particle n at the time t. The neighbors are 
determined through Voronoi tessellation [3] . The time depen¬ 
dence of the bond angle correlations was monitored by calcu¬ 
lating C^{t) = (E„ / (E„ l'l>?(0)n where * 

denotes the complex conjugate. 

To measure bond-orientational relaxation in 3D we de¬ 
fine Ql^{t) = where qim{9,(l)) 

are the spherical harmonics [23] and the sum is over the 
neighbors of a particle i at a time t determined through 
Voronoi tessellation. Next, we define the correlation func¬ 
tion Qiit) = ([(47r)/(2Z-hl)]EiEm=-/<3im(f)[Qim(0)]*). 
We calculated CQ{t) = Qe{t)/Qe{0) to monitor the decay 
of orientational correlations. 

We note that the conclusions remain unchanged if we de¬ 
fine neighbors as being less than a distance equal to the first 
minimum of the pair correlation function rather than through 
Voronoi tessellation. 
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Here we present results for one additional 3D glass former and three additional 2D glass formers to demonstrate 
that the main features discussed in the text is due to dimensionality and not due to a peculiar feature of a specihc 
system. The additional systems are discussed in Methods and we label the additional systems as follows: (1) the 50:50 
mixture with the e(cr/r)^^ potential we label 50:50, (2) the mixture with 32.167% small particles and the e[(j/rY‘^ 
potential we label 32:68, (3) the harmonic sphere system we label Harm. We also present some additional results for 
the 2D system discussed in the main text, the 65:35 Kob-Andersen mixture, and we label this system KA. Almost all 
the units are standard reduces units where (taa is the unit of length and caa is the unit of energy and the mass is 
one for all the systems. The only exception is that the temperature unit for the the harmonic sphere system is 10“^e. 
We label each hgure or curve in the hgure with the dimensionality. 






FIG. 1: Self intermediate scattering function and bond-orientational correlation function, (a)-(c) Fs{k;t) (solid 
lines) and Cq/{t) (dashed lines) for the 50:50 system (a), the 32:68 system (b), and the Harm system (c). (d) The ratio refra 
for the four 2D glass formers. 
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FIG. 2: Translational and bond-orientational relaxation comparison, (a)-(b) The self-intermediate scattering function 
Fs{k;t) (solid lines) and the bond-orientational correlation function rescaled by the alpha-relaxation time Tq, in 3D 

(a), and 2D (b). 




FIG. 3: Finite size effects, (a) The mean square displacement for the KA mixture in 2D for different system sizes at 
T = 0.45. (b) The self-intermediate scattering function for the standard Newtonian dynamics (solid lines) and the inherent 
structure dynamics (dashed lines) for the KA system in 2D for different system sizes at T = 0.45. 
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FIG. 4: Dynamic correlation length. The dynamic correlation length ^ 4 (tq) versus the alpha-relaxation time Tq, for the 
KA system (filled symbols) and the 32:68 system (open symbols) in 2D (green circles) and 3D (red squares). 



